function valentinesRose1
% 绘制玫瑰花
x=linspace(0,1,300);
theta=linspace(-2*pi,15*pi,300);
[x,theta]=meshgrid(x,theta);
phi=(pi/2).*exp(-theta./8./pi);
X=1-.5.*(1.25.*(1-mod(3.6.*theta,2*pi)./pi).^2-1/4).^2;
y=1.95653.*x.^2.*(1.27689.*x-1).^2.*sin(phi);
r=X.*(x.*sin(phi)+y.*cos(phi));
roseHdl=surf(r.*cos(theta),r.*sin(theta),X.*(x.*cos(phi)-y.*sin(phi)),'EdgeColor','none');

% 蓝白色
roseHdl.CData=r.*cos(theta);
CM=slanCM('blues');
colormap(CM)

% romao
% roseHdl.CData=r;
% CM=slanCM(134);
% colormap(CM)

% twilight
% roseHdl.CData=r;
% CM=slanCM(136);
% colormap(CM)

% 颜色列表插值配色
% roseHdl.CData=r;
% CM=[0.5300    0.8300    0.8100
%     0.5200    0.7500    0.8200
%     0.4900    0.6200    0.8400
%     0.4900    0.5600    0.8400
%     0.4700    0.4900    0.8500
%     0.4500    0.3500    0.8700
%     0.9500    0.9500    0.9500];
% CMX=linspace(0,1,size(CM,1));
% CMXX=linspace(0,1,256)';
% CM=[interp1(CMX,CM(:,1),CMXX,'pchip'),interp1(CMX,CM(:,2),CMXX,'pchip'),interp1(CMX,CM(:,3),CMXX,'pchip')];
% colormap(CM)


% 坐标区域修饰
ax=gca;hold on;grid on;
axis([-1,1,-1,1,-.5,1])
ax.FontName='Cambria';
ax.LineWidth=1;
ax.GridLineStyle='-.';
ax.Projection='perspective';
ax.XMinorTick='on';
ax.YMinorTick='on';
ax.ZMinorTick='on';
end